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A study of using an explicit nonlinear dynamic finite element code for simulating the 
water landing of a space capsule was performed. The finite element model contains 
Lagrangian shell elements for the space capsule and Eulerian solid elements for the water 
and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method 
were used for predicting the fluid and structure interaction forces. The space capsule was 
first assumed to be rigid, so the numerical results could be correlated with closed form 
solutions. The water and air meshes were continuously refined until the solution was 
converged. The converged maximum deceleration predicted is bounded by the classical von 
Karman and Wagner solutions and is considered to be an adequate solution. The refined 
water and air meshes were then used in the models for simulating the water landing of a 
capsule model that has a flexible bottom. For small pitch angle cases, the maximum 
deceleration from the flexible capsule model was found to be significantly greater than the 
maximum deceleration obtained from the corresponding rigid model. For large pitch angle 
cases, the difference between the maximum deceleration of the flexible model and that of its 
corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat 
shield qualitatively support the findings presented in this paper. 


I. Introduction 

Many practical engineering problems require detailed water landing impact analyses. Examples of such 
problems include ship slamming, torpedo water entry, and space capsule sea landing. Analytical approaches for a 
rigid body water landing developed by von Karman 1 and Wagner 2 in the early nineteenth century are still being 
widely used for obtaining closed form solutions. A review describing water entry studies related to aerospace 
structures between 1929 and 2003 was recently published 3 . This review pointed out that finite element techniques 
for analyzing both fluid and structural response in a single model were only recently available. It also stated that the 
aircraft industry requires a validated numerical tool in which the interaction between the fluid and the structure can 
be modeled accurately. One final suggestion of the paper is that numerical modeling techniques such as hybrid 
analysis, finite element modeling and smoothed particle hydrodynamics 4,5 will certainly be utilized for future work; 
however a large amount of work needs be done to fully validate these numerical modeling techniques. 

Recent advancement in the high-speed computing technology for large-scale analyses and the development of 
the advanced explicit nonlinear dynamic finite element method provides a powerful capability for solving the fluid 
and flexible structure interaction (hydroelastic) problems. Explicit nonlinear dynamic finite element codes with an 
Arbitrary Lagrangian Eulerian (ALE) solver 67 for modeling water landing impact are currently available in 
commercial codes such as LS-DYNA 8 8 and MSC.Dytran 69 . These codes have been used by both researchers and 
engineers for simulating water impact of aircraft, spacecraft, and ship hull structures 10 " 12 . These studies found that 
the computational methods implemented in the explicit nonlinear dynamic finite element codes can produce useful 
results and provide valuable insight into water impact problems. However, these studies also revealed the 
shortcomings of the codes such as the mesh dependency of the prediction results, the lack of good correlations 
between test data and numerical predictions, and the computationally intensive issues related to the explicit method 
and the ALE solver. 
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In the Apollo program, a vast amount of test data and analytical solutions related to water landing impact was 
developed and documented in many publications, such as Refs. 13 to 16. Water landing was used exclusively for 
returning the Apollo space capsule to the earth as shown in Fig. 1. The Apollo water landing with three parachutes 
deployed is shown in the figure. Most of the analytical solutions published in the Apollo reports were based on the 
von Karman approach, and the space capsule was considered to be a rigid body. However, some Apollo water 
landing test data showed that for certain impact orientations the bottom-structure’s flexibility had a significant effect 
on the maximum accelerations at the center of gravity (CG) of the space capsule. For these cases, the predictions 
based on the rigid body spacecraft often underestimated the actual structural responses 15,16 . 

The flexible bottom structure of a space capsule is the heat-shield that may be made of stiffened or sandwiched 
panels with non-uniform thicknesses. Therefore, closed form solutions for a flexible capsule water landing can be 
difficult to obtain. Since the explicit finite element codes are suitable for modeling complex flexible structures and 
fluid and structure interactions, it is expected that these codes will be extensively used in simulating the water 
landing of future space capsules. 

This paper provides the results of a preliminary study using LS-DYNA®, an explicit nonlinear dynamic finite 
element code, for simulating a conceptual space capsule water landing. The objective of this preliminary study is to 
demonstrate that the numerical results obtained are reasonable and can be used to adequately assess design 
configurations. Several modeling aspects were studied, including determination of an adequate mesh size; 
comparison of the rigid body results with classical solutions using the von Karman and Wagner approaches; and 
comparison of the numerical results of a flexible capsule model with published Apollo water impact test data. 

The following sections describe the finite element models, computational methods, and numerical results. A 
representative finite element model of a conceptual space capsule with the surrounding air and water is presented in 
Section II, the computational methods and approaches are given in Section III, a mesh refinement study is presented 
in Section IV, and the water landing simulations for both the rigid model and the flexible model are presented in 
Sections V and VI, respectively. At the end of the paper, concluding remarks are given in Section VII to summarize 
the findings from this study. 

II. Finite Element Model of a Conceptual Space Capsule 

An early conceptual design of an Apollo-like space capsule for future space exploration is shown in Fig. 2. This 
design has a maximum diameter of 198 inches and a height of 130 inches. The weight is estimated to be 16,200 lbs. 
The assumed geometries, CG locations, and mass moments of inertia are also shown in the figure. Two coordinate 
systems are shown in Figure 2: the one at the cone apex is the global coordinate system and the other at the CG is 
the local coordinate system. The local coordinate system is attached to and rotates with the vehicle. 

A representative finite element (FE) model is shown in Fig. 3. The capsule was discretized with 4-node quad 
shell elements. The nominal edge length of the capsule shell elements in the impact area is about 3.0 in. Coarser 
meshes were used for the other areas. The air and water were discretized with 8-node solid elements with an ALE 
multi-material element formulation 8 . In the interaction region shown in Fig. 3, the ALE solid element size is 2.4 in. 
x 2.4 in. xl.2 in. The air was modeled as a vacuum material. The water domain has 1,218,060 nodes and 1,180,000 
elements, and the air domain has 832,341 nodes, and 820,000 elements. The depths of the water and air are 256.8 in. 
and 168.0 in., respectively. 

Note that the mesh of the capsule is immersed in the Eulerian meshes of water and air, but the fluid nodes and 
structure nodes do not need to be coincident. Figure 3 shows the local coordinate system (fixed on the capsule) that 
was used for outputting all the CG velocity and acceleration results in this paper. These results can be compared 
directly with experimental results without coordinate transformations. Vy and Vh denote the positive directions of 
vertical and horizontal velocities, respectively. For all analysis cases performed in this study, the vertical landing 
velocity (Vy) is assumed to be 25 ft/s in the global X direction and the horizontal velocity (Vh) is assumed to be 0 
ft/s in the in the global negative Z direction. 

In this study, air was modeled as a vacuum and water was modeled as a null hydrodynamic material. An 
equation of state (EOS) was used to compute the water pressure 17 ls . For water in compression, the equation of state 
used is 
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(1) 


P(p)= K( — - 1) 

Po 

where P is the pressure, p is the instantaneous density, K is the initial bulk modulus, and p a is the initial density. 

The water is assumed to be slightly compressible during landing impact, so the high order terms in Reference 18 are 
neglected, and the water temperature remains constant. 

TIT. Computational Methods and Approaches 

The explicit finite element code, LS-DYNA", and an ALE solver were used to analyze the model shown in Fig. 
3, and a penalty method was used to determine the fluid-structure interaction forces (contact forces) between the 
capsule and the water. This penalty method, allowing fluid to flow around a structure but not through the structure, 
was used to determine the impact force on a fluid node based on the penetration depth of this node into the structure. 
The penalty forces were then applied to both the fluid and structure to push the fluid out of the structure. 

This study focused on two technical issues: mesh sensitivity and structural flexibility. Mesh refinement studies 
were performed first to determine an adequate water mesh size. A rigid capsule model was used for the mesh 
refinement study, so the results could be compared with the closed form solutions. Then, the water landing of a 
capsule with a flexible bottom structure, referred to as the flexible model, was simulated, and the results were 
compared with test data in the literature. A flow chart given in Fig. 4 depicts the approaches used for this study. 

IV. Mesh Refinement Studies 

To investigate whether a converged solution could be obtained through mesh refinements, FE models with a 
rigid capsule at a zero-degree pitch angle, using six different water mesh sizes (types) for the impact interaction 
region, were analyzed. Dimensions of the six different mesh types are listed in Table 1. The vertical landing velocity 
(V v ) is 25 ft/s, and the horizontal velocity (V H ) is 0 ft/s. 


Table 1 Solid element dimensions for various mesh types 


W L 

:>] D 

Solid Element 

Mesh Types 

L, in 

W, in 

D, in 

6x6x6 

6 

6 

6 

6x6x2. 4 

6 

6 

2.4 

2. 4x2. 4x2. 4 

2.4 

2.4 

2.4 

2 .4x2 .4x1. 2 

2.4 

2.4 

1.2 

1.2x1. 2x1. 2 

1.2 

1.2 

1.2 

1.2x1. 2x0.6 

1.2 

1.2 

0.6 


Results for the CG X-velocity and X-acceleration histories are shown in Figs. 5 and 6, respectively. For all the 
FE models, the bottom of the capsule was discretized with approximate 3 in. x 3 in. quadrilateral elements. In this 
paper, the negative value of acceleration is the deceleration; and the acceleration and the absolute value of 
acceleration are used interchangeably depending on the content of the text. Note that the predicted X-acceleration 
time histories contain high frequency oscillations. These high frequency oscillations were filtered out with a 100 Flz 
low-pass filter (FIR100) 19 . The filtered time histories are plotted in Fig. 6. Figure 5 shows that the CG velocity for 
Mesh Type 6x6x6 decreases much faster than that computed with the other finer mesh models. Comparison of the 
X-velocity curves for Mesh Types 2. 4x2. 4x1.2 and the 1.2x1. 2x1. 2 showed very small difference. Figure 6 shows 
that the maximum X-acceleration at the CG is continuously reduced as the mesh size becomes finer. Comparison of 
the X-acceleration curves for Mesh Types 2. 4x2. 4x1. 2 and 1.2x1. 2x1. 2 showed that nearly the same maximum 
accelerations are obtained. Note the results for Mesh Type 1.2xl.2x 0.6 are nearly the same as the results for Mesh 
Type 1.2x1. 2x1. 2 as shown in Figs. 5 and 6. For clarity, there are only a few data points for Mesh Type 1.2xl.2x 
0.6 shown in Fig. 5. 

Figure 7 contains the CG maximum X-acceleration as a function of the element volume. The mesh refinement, 
from a size of 2.4 in. x 2.4 in. x 1.2 in., to 1.2 in. x 1.2 in. x 1.2 in., and then to 1.2 in. x 1.2 in. x 0.6 in., yields no 
appreciable change in the maximum X-accelerations. A refinement of the capsule mesh was also performed in which 
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the quadrilateral elements at the bottom of the capsule were refined from approximate 3 in. x 3 in. to 2 in. x 2 in., 
and no noticeable change in results from the refined model was found as shown in Fig. 7. By examining the plot in 
Fig. 7, it is clear that a converged solution has been obtained. The analytical solutions 22 2 ' using von Karman and 
Wagner approaches are also shown in Fig. 7 to illustrate that the converged solution is bounded by the two 
analytical solutions. 


V. Analysis Results for Rigid Capsule Model 

Water landing analyses results for the capsule, assumed to be rigid, with various pitch angles are presented in 
this section. In these analyses, Mesh Type 2. 4x2. 4x1. 2 was chosen to model the fluid in the interaction region for 
computational efficiency. The water model is the same as the one shown in Fig. 3. The time step required for the 
water landing analysis is 6.84xl0" 6 s. A typical analysis time for simulating the first 50 ms of the water impact took 
about 28.5 CPU hours on a single 2 GFlz Opteron 1 ' 864 processor. 

To illustrate that the capsule position and the water splash (jet flow) can be quite realistically simulated, a snap 
shot at 0.15 s after impact of the capsule with a fifteen-degree pitch angle landing is shown in Fig. 8. The effect of 
jet flow on the impact force prediction was considered in Wagner’s solutions 2,20 . For this water landing case, the X- 
acceleration and the Y-rotational acceleration time histories are plotted in Figs. 9 and 10, respectively. Theses two 
types of accelerations can be combined for estimating the acceleration at other locations on the vehicle. The 
accelerations along the Y and Z axes are small and not presented in this paper. As shown in Figs. 9 and 10, the 
acceleration curves contain numerous high frequency oscillations. During the impact, the wetted surface of the 
bottom structure grows continuously with time, and the peak pressure is located at the edge of the wetted area 13,15 . 
In the explicit finite element simulations, the continuous growth of the wetted area introduces a series of discrete 
impacts, each of which is associated with a spike in the interface force 21 . These intermittent spike forces are 
considered to be the sources of the oscillatory noise in the acceleration history curves. The high frequency noise was 
filtered out with a FIR 1 00 11 filter. Filtered data are plotted in Figs. 9 and 10 for the X-acceleration and the Y- 
rotational acceleration, respectively. 

Filtered X-acceleration time histories for all pitch angles are plotted in Fig. 1 1 for comparison. For pitch angles 
less than or equal to 20 degrees, the maximum accelerations occurred within the first five milliseconds after impact. 
X-accelerations were small for pitch angles of 30 and 40 degrees. The maximum acceleration as a function of pitch 
angle is plotted in Fig. 12 which shows that the maximum accelerations for 0, 10, and 15-degree pitch angles are 
around 15 g and for 30 and 40-degree pitch angles are below 5 g. There is a large reduction of the maximum X- 
acceleration when the pitch angles increase from 15 degrees to 30 degrees. A similar trend was found in the Apollo 
reports reviewed. Note that an entry angle of 27.5 degrees was used for the Apollo space capsule water landing 15 . 
The analytical solutions 22, 2,1 using von Karman and Wagner approaches are also plotted in Fig. 12. Since the 
maximum X-acceleration for the 0-degree pitch angle is bounded by the two analytical solutions, this indicates that 
the numerical solutions obtained are reasonable. 

Y-rotational acceleration histories for all pitch angles are plotted in Fig. 13 for comparison. Since the CG is not 
located at the capsule’s axis of symmetry (see Figure 2 for the CG locations), small negative Y-rotational 
accelerations were predicated instead of zero for the 0-degree pitch angle. Large Y-rotational accelerations were 
found for pitch angles of 15 degrees and 20 degrees. For pitch angles above 30 degrees, the Y-rotational 
accelerations decreased to below 20 rad/s 2 . The maximum Y-rotational accelerations as a function of pitch angles 
are plotted in Fig. 14. The maximum Y-rotational acceleration curve reaches its peak about 60 rad/s 2 at a pitch angle 
around 17 degrees. 


VI. Analysis Results for Flexible Capsule Model 

For an actual capsule, the bottom structure that impacts the water is expected to be flexible. Experimental data 
presented in the Apollo reports 14 15 show that the interaction, between water and a flexible bottom structure (heat 
shield), can generate much larger impact forces than those of a rigid capsule resulting in larger CG acceleration 
levels. This increase in acceleration was found for low impact angles. For impact angles greater than 30 degrees, the 
effect of flexibility is negligible. This part of the study demonstrates that the hydroelastic effect between water and 
the capsule’s flexible bottom structure can be numerically simulated. 
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The flexible capsule model shown in Fig. 15 has a flexible bottom portion and a rigid upper portion. The flexible 
bottom portion, modeled using shell elements, is assumed to be a linear elastic material that has an equivalent 
Young’s modulus of 5.0xl0 5 6 psi, a thickness of 1.0 in., and a density of 0.1 lb/in’. The CG location, total mass, and 
moments of inertia of the flexible model are the same as the previous rigid body model. In the simulations, the 
capsule landed on water with a vertical velocity of 25 ft/s and at five pitch angles: 0, 10, 20, 30, and 40 degrees. This 
flexible model water landing analysis required a time step about one quarter of the time step required for the rigid 
body model water landing analysis. A typical analysis for simulating the first 50 ms of the water impact took about 
173.5 CPU hours on a single 2 GHz Opteron ! “ 864 processor. 

The interaction between the water and the capsule’s flexible bottom was investigated to understand the 
hydroelastic effect on the X-accelerations. For a 0-degree pitch angle landing case, the X-displacement histories for 
the vehicle CG and for the bottom center point are plotted in Fig. 16. At any instant, the difference between these 
two curves is the deflection at the center of the bottom surface. A maximum upward deflection (-X direction) of 
1.643 in. occurred at 14.75 ms. Note, at this instant, the CG and the bottom center point had the same velocity as 
shown in Fig. 17. After reaching the maximum deflection, the capsule’s bottom surface sprang back and created a 
downward velocity much higher than the initial capsule velocity. This maximum velocity is about 550 in/s at the 
bottom center point and is shown in Fig. 17. This high velocity interaction resulted in a large impact force which in 
turn caused a large deceleration (more than 30 g) of the capsule as shown in Fig. 18. Figure 19 shows that the 
velocity of the flexible model after 14.75 ms decreases at a much faster rate than that of the rigid model. 

The fdtered X-acceleration and Y-rotational acceleration histories of the capsule water landing with five 
different pitch angles are plotted in Figs. 20 and 21. These acceleration histories were compared with those for the 
rigid capsule shown in Figs. 11 and 13. The maximum acceleration for the flexible model was found to occur at a 
much later time, around 20 ms, for the same pitch angle landing. Figure 22 plots the maximum acceleration as a 
function of pitch angle for both the rigid model and the flexible model. The von Karman and Wagner solutions from 
the Appendix for the 0-degree pitch angles are also shown in the figure. It is apparent that the maximum 
accelerations for the flexible model at low pitch angles (< 20 degrees) are much greater than those predicted by the 
rigid model. However, the difference between the maximum accelerations from the flexible model and the 
maximum accelerations from the rigid model diminishes as the pitch angle increases. Maximum Y-rotational 
accelerations from both the rigid model and the flexible model are plotted in Fig. 23. For low pitch angles, much 
greater Y-rotational accelerations were predicted for the flexible model than for the rigid model. The maximum Y- 
rotational accelerations predicted for the flexible model at pitch angles greater than 20 degrees are nearly the same 
as those predicted with the rigid model. 


VII. Concluding Remarks 

A preliminary study of simulating the water landing of a conceptual space capsule with an explicit finite element 
code was presented. Tasks performed in this study include conducting a mesh refinement study, comparing the 
numerical solutions with analytical solutions for the capsule modeled as a rigid body, and investigating the 
acceleration magnified by the hydroelastic effect for a capsule with a flexible bottom structure. 

Mesh refinements were first performed, so an adequate mesh size could be identified for modeling the water and 
air for the rest of this study. Both the fluid mesh and the capsule mesh were refined to investigate whether a 
converged solution could be obtained. Since the change of the solutions was found to be insignificant when the 
meshes were sufficiently refined, this clearly indicates that a converged solution has been obtained. 

Analyses were performed for the rigid capsule model entering the water at various pitch angles. The maximum 
acceleration is about 15g for pitch angles less than 15 degrees and less than 5 g for pitch angles above 30 degrees. 
Also, it was shown that the maximum X-acceleration for the 0-degree pitch angle case is bounded by the closed 
form solutions based on the von Karman and Wagner approaches. This indicates that the analyses performed can 
produce satisfactory results to use in design studies. 

Water landing analyses of the flexible capsule model with various pitch angles were performed. The maximum 
X-accelerations for low pitch angles obtained from the flexible model are significantly greater than the rigid body 
solutions, nearly twice as large for the 0-degree pitch angle case. Furthermore, the fluid-structure interaction process 
during the impact was examined, and it was found that the bottom structure deformed upward at the initial impact. 
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Next, the center portion of the bottom surface moved downward to re-impact the water at a speed much greater than 
the initial impact velocity. The maximum X-acceleration occurred during this re-contact process. The maximum X- 
accelerations for pitch angles above 20 degrees predicted by both the flexible model and the rigid model are nearly 
the same. Thus, for large pitch angles, the water landing solutions from the rigid model may be sufficient for the 
conceptual design. These findings from the flexible bottom model qualitatively agree with the experimental data of 
Apollo capsules with a heat shield attached. 

In summary, this preliminary study demonstrates that the explicit nonlinear dynamic finite element method is 
suitable for simulating the water landing of space capsules. By using an adequate water mesh size and accurately 
modeling the flexibility of the bottom structure, satisfactory numerical results can be obtained. 

Appendix: Closed Form Solutions from von Karman and Wagner Approaches 

For a space capsule that has a spherical bottom and is assumed to be rigid, closed form solutions based on the 
von Karman 1 and Wagner 2 approaches are available for correlating with the results from the explicit finite element 
analyses. The von Karman approach is based on the momentum theorem and using an added virtual mass. The 
penetration depth is determined without considering the splash-up of water level. The Wagner approach uses 
rigorous fluid dynamic formulations and considers the effect of the splash-up water level on the impact force. 

Analytical solutions for a spherical bottom body impacting with water using the von Karmen method, presented 
in Ref. 22, are summarized here. The magnitude of the virtual mass for a spherical bottom body is 



%(2 R-b) 


(Al) 


where m is the virtual mass, p is the mass density of water, b is the water depth, and R is the radius of the spherical 
bottom. The instantaneous velocity, V, of the center of gravity of the rigid body is 

V = — = V 0 (\ + mg/W)~ 1 

dt (A2) 

where t is time after impact, V„ is the initial velocity, g is the gravitational constant, and W is the weight of the rigid 
body. By substituting Eq. Al into Eq. A2, the instantaneous velocity can be rewritten as 


v = 


V 0 

I , 8 pgR 2 , b 3/2 
3 W K R 


(A3) 


The overall acceleration, a, can be written as 


a = 


d 2 b 

dt 2 


3x2 


1/2 
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UJ 
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4npgR 


+ 2 /l x- 1 
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Assuming « 1 , the maximum acceleration can be found as 


(A4) 
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with the impact time at 


(A5) 
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and the penetration depth at 


(A6) 
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In the von Karman approach, the rise of water due to the splash up is not considered. The effect of splash up was 
considered by Wagner and found to have significant effect on the impact force. Recently, Miloh 2 ' used a semi- 
Wagner approach to determine the slamming coefficient, a non-dimensional parameter that is defined as 

C,{b/R)= 2F 


pxR% 2 


where F is the impact force. Base on his analytical derivations, Miloh proposed that 


(A8) 
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is suitable for initial stage slamming. The maximum acceleration can be estimated as 

* _ g 


2 W 


C s {b mwc /R)pxR 2 V 0 1 


(A9) 


(A10) 
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Fig. 1 Apollo water landing with three parachutes deployed. 
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Fig. 2 Dimensions of the space capsule, the global coordinate system at the apex, and a local coordinate 
system at CG. 



Fig. 3 Finite element model and the local coordinate system at CG with positive vertical velocity (V v ) and 
positive horizontal velocity (V H ) normal and parallel to the water surface, respectively. 
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Fig. 4 Flow chart showing the approaches used for this study. 



Fig. 5 CG X-velocity histories for various size meshes. 
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Fig. 6 Filtered X-acceleration histories for various size meshes. 



Fig. 7 Maximum X-acceleration converged as the water mesh size refined. 
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Fig. 8 Capsule water landing with 15-deg pitch angle showing the splash of water. 



Time, sec 

Fig. 9 Unliltered and filtered X-acceleration histories for the water landing of the rigid capsule with a 15- 
deg pitch angle. 
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Fig. 10 Unfiltered and filtered Y-rotational acceleration histories for the water landing of the rigid 
capsule with a 15-deg pitch angle. 



Fig. 11 Filtered X-acceleration histories for the water landing of the rigid capsule with various pitch angles. 
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Pitch Angle, deg 

Fig. 12 Maximum X-acceleration as a function of pitch angle for the water landing of the rigid capsule with 
various pitch angles. 



Fig. 13 Filtered Y-rotational acceleration histories for the water landing of the rigid capsule with various 
pitch angles. 
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Fig. 14 Maximum Y-rotational acceleration as a function of pitch angle for the water landing of the rigid 
capsule. 



Fig. 15 The flexible capsule model containing a flexible bottom portion and a rigid upper portion. 
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Time, sec 

Fig. 16 X-displacement histories at the CG and the center of the bottom structure for 0-deg pitch angle 
landing. 



Time, sec 

Fig. 17 Velocity histories at CG and the center of the bottom structure for 0-deg pitch angle landing. 
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Time, sec 

Fig. 18 Filtered and unfiltered X-acceleration histories for the water landing of the flexible model with 0- 
deg pitch angle. 



Fig. 19 Velocity histories for the rigid capsule and the flexible capsule models with a 0-deg pitch angle. 
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Fig. 20 Filtered X-acceleration histories for the flexible model with various pitch angles. 



Time, sec 

Fig. 21 Filtered Y-rotational acceleration histories for the water landing of the flexible model with 
various pitch angles. 
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Fig. 22 Maximum X-acceleration as a function of pitch angle for the water landing of the rigid and the 
flexible models. 



Fig. 23 Maximum Y-rotational acceleration as a function of pitch angle for the water landing of the rigid and 
the flexible models. 


19 

American Institute of Aeronautics and Astronautics 


